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Abstract 

We report a surprising hysteretic behavior in the dynamics of a simple one-dimensional 
nonlinear model inspired by the tribological problem of two sliding surfaces with a 
thin solid lubricant layer in between. In particular, we consider the frictional dy- 
namics of a harmonic chain confined between two rigid incommensurate substrates 
which slide with a fixed relative velocity. This system was previously found, by ex- 
plicit solution of the equations of motion, to possess plateaus in parameter space 
exhibiting a remarkable quantization of the chain center-of-mass velocity [dynamic 
pinning) solely determined by the interface incommensurability. Starting now from 
this quantized sliding state, in the underdamped regime of motion and in analogy to 
what ordinarily happens for static friction, the dynamics exhibits a large hysteresis 
under the action of an additional external driving force F ext - A critical threshold 
value F c of the adiabatically applied force F ext is required in order to alter the ro- 
bust dynamics of the plateau attractor. When the applied force is decreased and 
removed, the system can jump to intermediate sliding regimes (a sort of "dynamic" 
stick-slip motion) and eventually returns to the quantized sliding state at a much 
lower value of F eX f On the contrary no hysteretic behavior is observed as a function 
of the external driving velocity. 
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1 Introduction 



Nonlinear systems driven far from equilibrium exhibit a very rich variety of 
complex spatial and temporal behaviors [1]. In particular, in the emerging 
field of nanoscale science and technology, understanding the nonequilibrium 
dynamics of systems with many degrees of freedom which are pinned in some 
external potential, as is commonly the case in solid state physics, is becoming 
more and more often an issue. Friction belongs to this category too, because 
the microscopic asperities of the mating surfaces may interlock. It has been 
frequently shown [2] that simple phenomenological models of friction give good 
qualitative agreement with experimental results on nanoscale tribology or with 
more complex simulation data of sliding phenomena. In this kind of simplified 
approaches, studies are typically restricted to describing microscopic dynam- 
ics in one (ID) or two (2D) spatial dimensions. The substrates defining the 
moving interface are modelled in a simplified, although often effective, way as 
purely rigid surfaces or as one- or two-dimensional arrays of particles inter- 
acting through simple (e.g., harmonic) potentials. Despite this crude level of 
description, this approach has frequently revealed the ability of modelling the 
main features of the complex microscopic dynamics, ranging from regular to 
chaotic motion. 

One of the pervasive concepts of modern tribology - with a wide area of 
relevant practical applications as well as fundamental theoretical issues - is 
the idea of free sliding connected with incommensurability. When two crys- 
talline workpieces with lattices that are incommensurate (or commensurate 
but not perfectly aligned) are brought into contact, the minimal force re- 
quired to achieve sliding, i.e. the static friction, should vanish, provided the 
two substrates are stiff enough. In such a geometrical configuration, the lattice 
mismatch can prevent asperity interlocking and collective stick-slip motion of 
the interface atoms, with a consequent negligibly small frictional force. Experi- 
mental observation of this kind of superlubric and anisotropic regime of motion 
has recently been reported [3,4]. The remarkable conclusion of frictionless slid- 
ing can be drawn, in particular, in the context of the Frenkel-Kontorova (FK) 
model (see [5] and references therein). Since however the physical contact be- 
tween two solids is generally mediated by so-called "third bodies", the role 
of incommensurability has been recently extended [6] in the framework of a 
driven ID confined model inspired by the tribological problem of two sliding 
interfaces with a thin solid lubricant layer in between. The moving interface is 
thus characterized by three inherent length scales: the periods of the bottom 
and top substrates, and the period of the embedded solid lubricant structure. 
In particular, in the presence of a uniform external driving, the interplay be- 
tween these incommensurate length scales can give rise to intriguing dynamical 
phase locking phenomena and surprising velocity quantization effects [7,8]. 
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Extending a previous study of this confined tribological model [7], here we 
focus on the remarkable hysteretic behavior that this system exhibits starting 
now from this quantized sliding state. We find a strictly analogy to what ordi- 
narily happens for static friction in the underdamped regime of motion [9] . The 
lubricant center-of-mass (CM) velocity turns out to be robustly locked to the 
quantized plateau value (dynamic pinning) which is only abandoned above a 
critical force. As long as inertia effects are not negligible compared to dissi- 
pative forces, the adiabatic variation (increase and decrease) of the external 
applied force shows a large hysteresis loop in the V cm - F ext characteristics. 
Some differences between this dynamic locking and the usual static pinning 
are also briefly discussed. 



2 Confined model and numerical method 



Like in Ref. [7], we consider a simplified one-dimensional generalized FK model 
consisting of two rigid sinusoidal substrates, of spatial periodicity a + and a_, 
and a chain of harmonically interacting particles, of equilibrium length ao, 
mimicking the sandwiched lubricant layer, as schematically shown in the inset 
of Fig. 1. The two substrates move at a constant relative velocity V cxt = 
V- — V + . (In particular we set in full generality V + = and V_ = V cxt ). In 
order to probe the robustness of quantized dynamics of V cm , an additional 
constant force F ext is applied adiabatically to all chain particles. The equation 
of motion of the i-th lubricant particle becomes: 



1 

mxi = — - 



2tt ^ 2tt . . 
F + sin — Xi + F_ sin — (x, - V ext t) 



2 

+ K(x i+ i + Xi-i - 2xi) - -f(2xi - V cxt ) + F cxt , (1) 

where m is its mass. F± are the amplitudes of the forces due to the sinusoidal 
corrugation of the substrates. Presently we set F_/F + = 1 as the least biased 
choice. K is the chain spring constant defining the harmonic nearest-neighbor 
interparticle interaction. The penultimate damping term in Eq. (1) originates 
from two frictional contributions of the form —7 (xj — V+) — 7 (xj — VI), where 
7 is a viscous friction coefficient accounting phenomenologically for degrees 
of freedom inherent in the real physical system (such as substrate phonons, 
electronic excitations, etc.) which are not explicitly included in the model. 
The infinite chain size is managed - in the general incommensurate case - by 
means of periodic boundary conditions (PBC) and finite-size scaling (see, for 
example, Refs. [7,8]). We finally take a + = 1, m — 1, and F + = 1 as basic 
dimensionless units. 

The detailed behavior of the driven system in Eq. (1) depends crucially on the 
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relative (in)commensurability of the substrates and the chain. The relevant 
length ratios are defined by r± = a±/ao, we assume, without loss of generality, 
r_ > r + , focusing mostly on the case r + > 1. In particular, in order to make 
a comparison with previous tribological studies [7,8], we shall restrict our 
present analysis to the incommensurate golden-mean case = (v^ + l)/2« 
1.6180 with ratios (r+, r_) = (0, 2 ). Since the qualitative features of velocity 
quantization phenomena were proved [7] to survive for much more general 
values of r + and r_, this specific choice of incommensurability should not be 
considered too restrictive. 

The equations of motion (1) are integrated using a standard fourth-order 
Runge-Kutta algorithm. The system is initialized with the chain particles 
placed at rest at uniform separation a . After relaxing the starting config- 
uration and selecting a reference frame in which the bottom substrate is at 
rest (V + = 0), the top substrate starts sliding at the imposed constant veloc- 
ity V- = V cxt . For a very wide range of model parameters the system reaches, 
after an initial transient, the quantized "dynamical stationary" state with 
Km/Koct = V'piatcau/Kxt = 1 — r^ 1 , that depend solely on the chosen incom- 
mensurability ratio r+. In order to investigate the possibility for the system 
to exhibit hysteresis starting now from this quantized sliding state, at which 
the confined layer is robustly pinned, the additional constant force F ext acting 
on all chain particles is now varied upward and downward adiabatically. 



3 Results and discussion 

After stationarity has been reached, Fig. 1 shows the striking behavior of 
the normalized time-averaged CM velocity of the sandwiched lubricant chain, 
Vcm/Vext, as a function of the stiffness K, for the two incommensurate case of 
golden mean (GM) and spiral mean (SM) 1 . 

As discussed in previous works [7,8], the crucial feature is the presence of 
perfectly flat V cm /V cxt plateaus, whose precise value (1 — r^ 1 ) is independent 
not only of K, but also of 7, V cxt , and even of F_/F + . Their occurrence 
was ascribed to the intrinsic topological nature of this quantized dynamics. 
The phenomenon is explained by one confining substrate rigidly dragging the 
topological solitons (kinks) that the embedded chain forms with the other 
substrate. 

Let us now turn to consider specifically the GM case only. Fixing a value of 

1 The golden mean = (s/b + l)/2 is the solution of the quadratic equation 2 — 
— 1 = 0; the spiral mean, a f» 1.3247 satisfying the equation a 3 — a — 1 = 0, 
belongs to the class of cubic irrationals. 
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the chain stiffness K lying approximatively in the middle of the plateau of 
Fig. 1 and considering a sufficiently small value of the damping coefficient 7 
(underdamped regime), we start investigating the hysteresis by applying an 
external force F ext to all the chain particles through an adiabatic increase and 
decrease process. 

The results are displayed in Fig. 2 for two different external driving veloci- 
ties. A clear hysteretic loop emerges, with qualitative similar features for high 
(upper panel) and low (lower panel) V cxt . Surprisingly, the cycle is broader for 
larger velocities. We will return to this point later on. 

The finding of exact plateaus implies a kind of "dynamical incompressibility" , 
namely, identically null response to perturbations or fluctuations trying to de- 
flect the CM velocity away from its quantized value. What is now the effect of 
the additional force F ext ? We find that as long as F cxt remains below a criti- 
cal threshold F c , it does perturb the single-particle motions but has no effect 
whatsoever on V cm , which remains exactly pinned to the quantized value, as 
could indeed be expected of an incompressible state. This picture is analogous 
to the pinning-depinning transition in static friction, where a minimum force 
(the static friction) is required in order to start the motion. Thus the sudden 
change of V cm taking place at F ext = F c can be termed a "dynamical depin- 
ning" . The value of F c is a nontrivial function of the parameters, and vanishes 
linearly when K approaches from below the upper border K c of the plateau. 
The depinning transition line F c , ending at K = K c , appears as a "first-order" 
line, with a jump AV^ in the average V cm and a clear hysteretic behavior as F 
crosses F c . As expected, we find that AV decreases to as K increases towards 
K c (not shown). Thus K = K c represents a sort of non-equilibrium critical 
point, where the sliding chain enters or leaves a dynamical orbit. The precise 
value of K c depends on parameters such as V cxt and 7; however, its properties 
do not. Fig. 3(a) displays the general decreasing behavior of K C) and thus the 
corresponding diminishing extension of the plateau, as a function of the exter- 
nal driving V cxt (no applied force F cxt ). As shown in Fig. 3(b), and contrary to 
what could be intuitively expected, no straightforward relation seems to exist 
between the plateau extension in K and its robustness against the external 
perturbing force F ext . At a fixed chain stiffness, quite large values of F c are 
found even for very high sliding velocities, where the plateau extension has 
already been reduced significantly. 

Depending on model parameters such as chain stiffness and external driving 
velocity, the dynamical depinning off the quantized sliding state takes place 
through different kinds of mechanisms, ranging from a series of intermittencies 
with a well-defined temporal periodicity to more chaotic and irregular jumps 
[7]. As displayed in Fig. 4, when the applied force is decreased and removed, 
the system may jump to intermediate sliding regimes and eventually returns 
to the quantized sliding state at a much lower value of F ext . As for the forward 
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dynamical depinning transition, the details of these backward steps in the 
Vcm ~ ^ext characteristics strongly depend on the parameters of the model. 
At high values of -F ext (= 0.07, left panels) the confined layer moves almost 
freely with a large sliding velocity depending on the damping coefficient 7. 
For intermediate force (F ext = 0.045, middle panels), simulations reveal the 
intriguing occurrence of a sliding regime closely reminiscent of a dynamic 
stick-slip motion. This intermittent dynamics is seen with particular clarity 
by plotting the particle trajectories in the reference frame which slides at the 
quantized velocity value of the plateau. A further reduction of F ext (= 0.02, 
right panels) brings the system back to the time-periodic dynamics of the 
quantized sliding state. 

The above picture of dynamical depinning as a first order transition is valid 
for weak dissipation. For strong dissipation, when the viscous damping coef- 
ficient 7 is much larger than the characteristic vibrational frequencies of the 
system (overdamped motion) the dynamical depinning is likely to be of sec- 
ond order: The forward and backward trajectories become indistinguishable, 
and hysteresis disappears, as shown in Fig. 5(a). In this strongly dissipative 
regime, we found instead of the hysteretic jumps a nonlinear mobility region 
of V cia versus F cxt , but without any bistability phenomenon. 

Finally, independently of the value of the chain stiffness, no hysteresis has been 
found in the underdamped regime by varying the external driving velocity with 
time with a gentle enough rate of increase and decrease. Panels (b) and (c) 
of Fig. 5 show, for two different K inside the quantized plateau region, the 
nonlinear, but not bistable, behavior of V cm /V cy ± as a function of adiabatic 
variation of V cxt . 



4 Conclusions 

We have shown that starting from the quantized V cm sliding state, previously 
found for a simple tribological model of a confined layer, the layer sliding dy- 
namics exhibits a large hysteresis under the action of an additional external 
driving force F ext trying to change V cm away from its quantized value. In anal- 
ogy to depinning in ordinary static friction, the hysteretic dynamical behavior 
depends strongly on whether the system degrees of freedom have sufficient 
inertia (underdamped regime) or if, on the contrary, the inertia is negligible 
(overdamped regime). 

The robustness of quantized dynamics is proved by the existence of a finite 
critical threshold F c needed to move the chain CM velocity away from the 
plateau value (dynamical depinning). When the applied force is decreased 
and removed, the system may jump to intermediate sliding regimes (a sort 
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of dynamic stick-slip motion) and eventually returns to the quantized sliding 
state at a much lower value of F cxt . 

There are however nontrivial differences from static friction. The first is that 
the dynamical pinning hysteresis cycle may be larger in situations where the 
pinning itself could be intuitively considered more fragile, e.g., for larger ex- 
ternal velocity. Another feature (presently under investigation, not discussed 
above) is that the sudden application of an external force can leave V cm locked 
to the quantized value, even if the applied force is larger than the dynamic 
depinning threshold F c , obtained instead through the adiabatic procedure 
sketched above. Once again, this is different from static depinning, usually 
requiring smaller force (than F s ) if applied suddenly. 

A final open question concerns the effect of a finite temperature on the dy- 
namical hysteresis. Preliminary results obtained through a Langevin dynamics 
indicate that, so long as the thermal energy is much smaller than an effective 
dynamical barrier (gap) preserving the incompressible plateau state, the ve- 
locity quantization is still observed. We expect, therefore, that the qualitative 
dynamical hysteretic behavior should not change too much. These aspects are 
currently under investigation. 
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Fig. 1. (Color online) Normalized velocity of the center of mass, V cm /V ex t, as a 
function of the chain stiffness K, for the golden mean (r + ,r_) = (0, (f> 2 ) and spiral 
mean (r + ,r_) = (a, a 2 ) incommensurability. Here 7 = 0.1 and V ex t = 0.1. Note 
the logarithmic scale in the abscissa. A sketch of the driven 3-length scale confined 
model is shown in the inset. 
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Fig. 2. (Color online) Hysteresis in the V cm — F cxt characteristics for the GM case and 
a relatively soft {K = 4) confined chain. The behavior is shown for high (V eX t = 0.1, 
upper panel) and low (V ex t = 0.01, lower panel) applied driving velocities. Adiabatic 
increase and decrease of F ext is denoted by triangles and circles, respectively. A 
characteristic multi-step feature appears when decreasing adiabatically the external 
force. 
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Fig. 3. (Color online) (a) (V ex t', K) phase diagram, with no additional force applied, 
for the GM incommensurability. At a fixed value of the external driving velocity, K c 
marks the point above which the quantized sliding regime of the chain breaks down, 
(b) Irregular dependence of the dynamical depinning force F c upon the driving 
velocity V ext for chain stiffness K = 4. 
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Fig. 4. (Color online) Time evolution of the velocities V{ (upper panels) and of the 
corresponding rescaled coordinates Xj— (V^iateau^) (lower panels) of seven contiguous 
chain particles. The plots refer to three different dynamical regimes observed in 
the adiabatic decrease process of the external force: free sliding at -F ext = 0.07 
(left), dynamic stick-slip at F ext = 0.045 (middle), and quantized sliding state at 
-Fcxt = 0.02 (right). Note the different scale in the ordinate of the t^-plots. Here 
7 = 0.1, K = 4, and V ext = 0.1. 
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Fig. 5. (Color online) No hysteretic behavior of V cm /V ex t is observed in leaving 
the V cm quantized plateau as a function of: (a) the external applied force in the 
overdamped regime; (b)-(c) the driving velocity, for different values of K. 
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